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O Abstract 

' We consider the problem of switching off unwanted interactions in a given multi- 

partite Hamiltonian. This is known to be an important primitive in quantum infor- 
mation processing and several schemes have been presented in the literature to achieve 
this task. A method to construct decoupling schemes for quantum systems of pairwise 
interacting qubits was introduced by M. Stollsteimer and G. Mahler and is based on 
orthogonal arrays. Another approach based on triples of Hadamard matrices that are 
closed under pointwise multiplication was proposed by D. Leung. In this paper, we 
show that both methods lead to the same class of decoupling schemes. Moreover, we 
establish a characterization of orthogonal arrays by showing that they are equivalent to 
decoupling schemes which allow a refinement into equidistant time-slots. Furthermore, 
we show that decoupling schemes for networks of higher-dimensional quantum systems 
with i-local Hamiltonians can be constructed from classical error-correcting codes. 



1 Introduction 

An important task in the study of quantum systems is to manipulate a given system Hamil- 
tonian by applying external control operations in such a way that in effect the time-evolution 
of some other desired target Hamiltonian is simulated. Typically, the available control op- 
erations are restricted and furthermore the control schemes employed are required to be 
efficient. Hence, the number of control operations should be a polynomial function in the 
number of particles which are governed the system Hamiltonian. In the context of pair- 
interaction (also called two-local or two-body) Hamiltonians acting on n qudits a repertoire 
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of techniques has been developed to use any entangling Hamiltonian for universal simulation 
of arbitrary couplings PQEHBHUISlinilZlIB]- Here the external control operations available 
are given by strong pulses which are local unitary rotations applied to the individual nodes. 
A cornerstone of this theory is the development of decoupling schemes and selective cou- 
pling schemes. Both are pulse sequences that switch off all unwanted interactions in a given 
Hamiltonian. In the case of a decoupling scheme all interactions have to be switched off. 
In contrast the requirement for a selective coupling scheme is that all interactions except 
for the interaction between two fixed nodes have to be switched off. Two methods have 
been proposed to achieve decoupling and selective coupling of a general pair-interaction 
Hamiltonians in quantum systems consisting of n qubits: 

Construction I This method, which was proposed by D. Leung [3], uses triples S x , S y , 
S z of Hadamard matrices. If the rows of these matrices satisfy a compatibility condition, a 
sequence of pulses around the a x , a y , and a z -&x.es in the Bloch sphere can be constructed. 

Construction II In the method put forward by M. Stollsteimer and G. Mahler [2j the 
pulses are constructed by using orthogonal arrays which are matrices which fulfill a bal- 
ancedness conditions between the rows. 

The purpose of this paper is to show that these constructions are equivalent, i.e., that 
each admissible triple of Hadamard matrices used in Construction I leads to an orthogonal 
array which can be used for Construction II (and vice versa). We first show this correspon- 
dence for systems consisting of n two-dimensional systems and generalize this in the sequel 
to higher-dimensional systems. Also the requirement that the given Hamiltonian to be a 
two-body Hamiltonian can be relaxed: we show that orthogonal arrays of strength t can be 
used to decouple any t-local Hamiltonian. 

2 The Framework: Average Hamiltonian Theory 

Switching off unwanted interactions is an important primitive in the approaches to render 
a given Hamiltonian to simulate any other Hamiltonian p] 121 E3 IU 03 El ■ Here simulation 
is usually understood in a narrow sense in which the desired target Hamiltonian is approx- 
imated up to terms of quadratic and higher orders. In the following we briefly introduce 
the facts of this framework of average Hamiltonian theory 111] which will be needed to 
develop the theory of decoupling schemes. 

Assume that the system Hamiltonian acts on an ra-fold tensor product Hilbert space 
H := C d ® C d ® . . . ® C d , where each C d denotes the Hilbert space of a so-called qudit. 
Let B := {a a \ a = 1, . . . , d 2 — 1} be a basis of traceless matrices acting on C d . The most 
general i-local Hamiltonian for a system of n coupled qudits is given by 

t d 2 -i 

S=l (kl,...,k s ) «l,-i«s = l 
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where the second sum runs over all s-tuples with (different) entries from {1,... , n} and 
Jk 1 ,...,k s -,a 1 ,...,a s £ C. Here and in the following we use to denote the operator that acts 
as A on the fcth qubit, i. e., A^ :=l®---®l®j4.<8)l<g)---®l. 

In the setting discussed in this paper the only possibilities of external control are given 
by local unitaries on each qudit. We assume that it is possible to implement them in- 
dependently. Formally, all control operations are elements of some finite subset C of the 
group I4(d)® n , where 14(d) denotes the group of unitary matrices acting on C d . A useful 
approximation is to assume that all operations in C can be implemented arbitrarily fast 
("fast control limit"). The simulation of Hamiltonians is based on the following average 
Hamiltonian ^01 approach. Let ti,t2, ■ ■ ■ ,t n be real numbers and V\, V2, ■ ■ ■ , Vjy £ C be 
control operations. Note that letting the system evolve for some time t (in the following 
referred to as "wait") has the effect to apply the unitary operator exp(-iHt). Hence the 
sequence 

perform V\, wait t\, perform V2, wait ... perform Vjy, wait ijv 

implements the evolution Yif=i ex P( — UjHUj tj), where Uj = Yi{=i^j- We say that the 
scheme consists of N intervals — sometimes also referred to as time-slots — and use the 
shorthand notation (t\, U%; t2, U2; • ■ ■ ; tjv> Un), where we tacitly assume that the underlying 
Hamiltonian H is fixed. If the times tj are small compared to the time scale of the natural 
evolution according to H this gives an approximation to the average Hamiltonian 

N 

where r := ^ ■ tj is the slow down factor, i. e., the relative running time of the evolution. 

We next introduce decoupling schemes which can be used to simulate the zero Hamil- 
tonian. For this reason they are also used in dynamical suppression of decoherence in open 
quantum systems ("bang-bang" control), see [12[ ll3 ( IT3j. A decoupling scheme is a sequence 
of control operations such that the resulting average Hamiltonian is the zero matrix for all 
system Hamiltonians of the form in eq. I|T|). 

Recall that a unitary operator basis (also called unitary error basis ^3]) is a collection of 
d 2 unitaries f7, that are orthogonal with respect to the inner product (A\B) := l/dtv(A^B). 
Bases of unitaries with this property were already studied by Schwinger \16^ and recently 
several explicit constructions have been found ^3 El DB| ■ 

Definition 1 (Decoupling scheme) A decoupling scheme D := (pi,U\; . . . ;pn,Un) is 
given by positive real numbers pj summing up to 1 and control operations Uj £ C such that 

N 

J2PjU}HUj = (2) 

j=i 

for all t-body Hamiltonians acting on n qudits. We call a decoupling scheme D regular if 
the lengths of the time-slots are the same, i.e., if pi = P2 = ■ ■ ■ = Pn an d in addition if the 
operators applied to each node form a unitary error basis. 
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Note that if the system consists of one d-dimensional node only, a decoupling scheme 
is equivalent to a unitary operator basis. This definition includes decoupling schemes con- 
sisting of time-slots of different length. Many of the decoupling schemes considered in the 
literature use only time-slots of equal length EB EJOi]- However, also decoupling schemes 
are used in which the intervals have different lengths, most notably the famous WaHuHa 
sequence |T3I ITO] . 

The difficulty for systems consisting of more than one node is that we still want to use a 
number of operations which is polynomial in the dimension d of the individual nodes as well 
as in the number n of nodes. To construct schemes with this property it is necessary to be 
able to apply selective pulses to the nodes [2] ■ In the following we present two constructions 
of schemes which achieve decoupling for any pair of nodes, i. e, these schemes can be used 
to decouple any pair-interaction Hamiltonian. By applying the same sequence of pulses to 
a fixed pair of nodes all interactions will be switched off with the exception for the bipartite 
system consisting of these two nodes. This in turn can be used for universal simulation. 



3 Hadamard Matrices, Sign Matrices, and Phase Matrices 

In the following we give a short account of the combinatorial objects underlying the con- 
struction used in [I] to obtain decoupling schemes for pair-interaction Hamiltonians acting 
on n qubits. The construction relies on the concept of so-called sign matrices which gener- 
alize the refocusing schemes for spin echo experiments on n qubits. The latter have been 
proposed in [Q and are based on Hadamard matrices. Since we will need Hadamard matrices 
for the subsequent constructions we briefly recall their definition. 

Hadamard matrices We denote the transposed of a matrix A by A t . A Hadamard 
matrix of order TV is a ±1 matrix Hjy of size N x N with the property that HnH 1 n = 
-/Vljy. Hadamard matrices have been studied in combinatorics for a long time and several 
constructions have been found. We refer to |21[ l22 | l23] for background on and constructions 
of Hadamard matrices. We give some examples for Hadamard matrices of small order (here 
and in the following the entries ±1 have been abbreviated to +/—): 



H 2 = , H 2 ® if 2 



It is known that a necessary condition for the existence of a Hadamard matrix is that either 
N = 2 or N = mod 4. A long-standing conjecture is whether indeed for any N = mod 4 
a Hadamard matrix of order N exists (2^- Since H 2 n ■= H 2 ® ... <8> H 2 (n tensor factors) 
is a Hadamard matrix, we obtain that in dimension N = 2 n at least one Hadamard matrix 
exists. 
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Sign matrices A sign matrix S n n of size n x N is given by the first n rows of a Hadamard 
matrix of order N. Hence S n! N is a ±1 matrix which satisfies S n ^S^ N = Nl n . Recall 
that the Schur product of two n x N matrices A and B is denoted by C := A o B and is 
defined as the entry-wise product: Cij := AijBij. As an example we define the following 
three sign matrices S x , S y , S z of size 7x8: 



Sx '■- 



/+-+-+-+-\ 
++--++-- 
+--++--+ 

+ + + + 

+-+--+-+ 
+ + + + 

V+--+-++-7 



Si I 



/+-+--+-+\ 
+-+-+-+- 

+ + + + 

++--++-- 
+--+-++- 
+-+--+-+ 

V+ + + +J 



, s z 



(+ + + + \ 

+--++--+ 
+--+-++- 

+ + + + 

++--++-- 
+-+--+-+ 

K+-+-+-+-J 



Besides the fact that they are closed under Schur product, i. e., S x o S y = S z , these matrices 
have another remarkable feature: all of their rows are actually rows of H2 (8> H2 <S> B<i- 
This guarantees that S x , S y , and S z are sign matrices of size 7x8. As we shall see in 
the next section when we study criteria for decoupling, these matrices cannot be used to 
decouple a general pair-interaction Hamiltonian on seven qubits since all rows of all three 
matrices together are not orthogonal (indeed, any row appears in each of the three matrices). 
However, if the Hamiltonian of a seven qubit network is of the particular form where only 
Ox <8) 0$ , o y k ^ ® Oy , and a z k ^ ® 0$ interaction terms occur in eq. (J2) , these matrices can 
be used for decoupling and selective coupling. 

In Section we will present criteria for decoupling pair-interaction Hamiltonians acting 
on qubits and show that they are fulfilled in case we can find sign matrices S x , S y , and 
S z which are closed under taking the Schur product and have the additional property that 
their rows are pairwise orthogonal. The approach 0j requires such triples S x , S y , and S z of 
sign matrices of size n x N which are related by S x o S y = S z . 



Phase matrices The restriction to consider orthogonal matrices with entries ±1 can be 
relaxed by allowing the entries to be more general complex phases. This gives additional 
flexibility for the decoupling of pair-interaction Hamiltonians acting on higher dimensional 
systems (qudits) and leads to the concept of phase matrices which are defined as follows: 
Let k £ N and let to = exp(2iri/k) G C be a primitive kth. root of unity. Then a phase matrix 
Pn t N of order k is an n x N matrix with entries in {l,co, ■ ■ ■ , u> k ~ 1 } such that P n ,NPn m = 
Nl n . Like in case of sign matrices, we are interested in collections of phase matrices 
which satisfy certain compatibility conditions. These conditions can be conveniently stated 
in terms of characters of some finite abelian group G. For the necessary background on 
characters of abelian groups see AppendixEI In the following we assume that the elements 
of G are given in a fixed order gi, ... , giQi and that the irreducible characters of G are in one- 
to-one correspondence with the elements of G and are given by {x g '■ g £ G}, cf. Theorem 
1151 Recall that the exponent e(G) of G is the smallest positive integer such that g e ^ = 1 
for all g G G. Now, let Pi, ■ ■ ■ ,P\g\ be phase matrices of order e(G) which are labeled by 
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the elements of G. We say that the P g , where g E G, are compatible with respect to the 
Schur product if 

PgOP h = P gh (3) 

holds for all pairs g,h £ G. Note that according to this definition sign matrices are a special 
case of phase matrices. Indeed, we obtain that any sign matrix is a phase matrix for the 
group G = 7Li x 7Li where the phase matrix corresponding to the identity is given by the 
all-ones matrix. 

The connection to decoupling schemes is as follows: the phase matrix Ph will describe 
the phase factors which are acquired when a fixed unitary matrix Uh is conjugated by some 
other matrices. The condition in eq. (j3J) ensures that each vector v^j := [Ph;k,j]heG f° r 
fixed k = 1, . . . , n and j = 1, . . . , N is a homomorphism from G to C x . Since the characters 
form a group and the order of the elements is fixed we obtain that v^j is a row from the 
character table of G, i. e., Vf-j = Xg f° r an element g £ G. Each element g £ G corresponds 
to a control operation U g . These facts will be used in Section [4.21 to show the vanishing of 
terms in a pair-interaction Hamiltonian in case a decoupling scheme is applied. 

For general G it turns out to be a non-trivial task to construct |G|-tuples of phase ma- 
trices which at the same time fulfill condition (f3j). In Section f5. II we will give a construction 
which is possible in case the dimension of the nodes is a prime power, i. e., |G| = p m , where 
p is prime and m £ N. 



4 Constructing Decoupling Schemes 

We continue the investigation of decoupling schemes with an observation concerning the 
relative lengths of the time-slots in the scheme. Recall that according to Definition ^ a 
scheme is regular if all intervals are of equal lengths. In Section 14.11 we show that not all 
decoupling schemes are regular, and that this is the case even if we are allowed to reorder and 
refine the time intervals. In Sections 14. 2l and l4.3l we will then introduce the two constructions 
for schemes for decoupling and selective coupling mentioned in the introduction. 



4.1 Decoupling schemes which are not regular 



Assume that a decoupling scheme D = (pj, Ui, P2, U2, ■ ■ ■ ,Pn-> Un) on a system consisting 
of n nodes is given. Since Uj is a local operation for each j = 1,...,N we have that 

Focusing on the first two nodes only, we can always obtain a new 



U; 



(«) 



decoupling scheme, which has the form 









sU {2) 
su N 





sU[ 2) 




sU {2) 
su N 



ru N 






su N 



where r, s £ M are such that rN = sN 2 



1. However, in order to obtain a 



regular scheme in general we cannot continue in this fashion to more than three nodes. The 
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Figure 1: A decoupling scheme for a system of three qubits which is not regular, i. e., the 
time-slots cannot be rearranged into a form where all time-slots have the same length. The 
transformations applied to the individual qubits correspond to the Pauli matrices as follows: 
1 = 1 2 , 2 = a x , 3 = a y , and 4 = a z . The time intervals indicated in the figure have four 
different basic lengths t\ , t% = 1/Ati, t 3 = -^t 2 , and £4 = v ^ 1 t 2 . For instance, the first 
intervals applied to the first qubit all length t±, whereas the interval lengths for qubit three 
are given by t 2 , t 2 , h, t 4 , t 4 , t 3 , etc. 



pulse sequence given in Figure H provides an example of a decoupling scheme that cannot 
be refined into time-slots which have the same lengths. The control operations used in the 
scheme are the Pauli matrices, which form a basis for the vector space of all 2 x 2 matrices 
and are given by 

1 \ / -i \ ( 1 









i 







1 j ' " y \ % y °* v - 1 

In the given example the network consists of three qubits and in each time-slot precisely 
one of the four Pauli matrices is applied. Indeed, we first verify that the pulse sequence 
defined in Figure ^ defines a decoupling scheme for any pair-interaction Hamiltonian: first, 
note that the sum of the times for each Pauli operator applied to the individual qubits 
is constant, i.e., the local terms are removed. Moreover, by considering pairs of rows we 
verify directly that also any pair of symbols (a, b) with a, b G {1,2,3,4} is applied for the 
same time t 2 . For example in case of rows two and three we obtain for the pairs (1,1) the 
total time t 2 and for (3,4) the total time £3 + t± = t 2 . However, the sequence cannot be 
subdivided into a finite number of intervals of equal lengths. Indeed, this would contradict 
the fact that y/2 is not a rational number. 



4.2 Decoupling schemes from sign and phase matrices 

We now describe the approach of |I] to construct decoupling and selective coupling schemes 
for general pair-interaction Hamiltonians acting on qubits. Then we present a generalization 
for qudits by generalizing the underlying group-theoretical structures based on our definition 
of phase matrices. 
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The qubit case: decoupling schemes from sign matrices A general pair-interaction 
Hamiltonian, i.e., two-local Hamiltonian, for n qubits may be written in the form 

h ■■= E E M^af + EE J ^ k) ( 4 ) 

k<£ a/3 k a 

where Jki;af3 £ Jk;a £ K and where a a are the Pauli matrices, i.e., a £ {x,y,z}. To 
construct decoupling schemes we choose I2, cr x , a y and cr 2 as control operations. Then 
in each time-slot of a decoupling scheme each of the n qubits is conjugated by precisely 
one of the matrices I2, cr x , Cy, and a z . Since the terms of the Hamiltonian in eq. (0J) are 
expressed in terms of the Pauli matrices we have to compute the resulting effect by this local 

(k) (k) (k) (k) 

conjugation: the possible sign assignments for lg , cr x , a y , a\ are given by the following 
table: 





1 2 




Oy 




1 2 
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+ 
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+ 
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+ 






Oy 


+ 










+ 






+ 



Here the rows are labeled by the operators used in eq. © and the columns by the operators 
realizing the conjugation at each time-slot. Note that in each column the signs multiply 
to +1 since conjugation by a local unitary corresponds to a SO(3) rotation of the Bloch 

(k) 

vector. Hence, the signs assignments of the three Pauli matrices o a acting on the same 
qubit k are not independent but rather the third is always given by the other two. 

(k) 

Now, in each interval a a acquires either a + or a — sign, which is controlled by the 
applied local unitaries (the identity matrix or a Pauli matrix). According to table © the 
bilinear coupling Jke;a/30a *crg is unchanged (negated) when the signs of or o4 agree 
(disagree) . 



Decoupling criteria in terms of sign matrices We show that sign matrices satisfying 
certain orthogonality conditions yield decoupling schemes. The (k,j) entry of S a for a G 
{x, y, z} is denoted by S a; kj and gives the sign of a a k ^ in the jth. time-slot. Hence a regular 
decoupling scheme which uses N time-slots can be obtained from these matrices using the 

following rules: if the triple of entries at position (k,j) is given by (+ + +), (H ),( — I — ), 

respectively ( h) then the operation applied to qubit k in time step j is given by I2, 

a x , <jy, respectively a z . 

Decoupling is achieved if any two rows taken from S x , S y , S z are orthogonal. To achieve 
selective coupling between two nodes k and £, the operations applied to the nodes k and t 
are chosen to be identical while still maintaining orthogonality of the modified sign matrices 
jlj. Within this framework sufficient and necessary conditions for decoupling of all pair- 
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interactions are given by the following equations: 



N 



(6) 



for all a and all k, and 



N 



S a; kjSp ; £j — (7) 

i=i 

for all a, (3 and all k < t. The first condition ensures that all local terms are removed and 
the second condition that all bilinear terms are removed. These conditions are satisfied if 
the sign matrices S x , S y , S z and all rows of all three matrices are orthogonal to each other. 

Example 2 As an example consider the following sign matrices which specify a decoupling 
scheme for a system of 5 qubits with 16 time-slots. 
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(8) 



(9) 



(10) 



In general the construction of admissible triples (S x , S y ,S z ) of sign matrices proves to be a 
delicate task. However, a construction of triples of sign matrices of size (2 2n — l)/3 x 2 2n by 
partitioning the rows of the Hadamard matrices Hf 2n , where n E N, was given in [3]. This 
construction has been revisited in where an alternative proof based on spreads in a finite 
geometry has been given. We will give yet another proof of this family of sign matrices in 
Section f5. II which is based on Hamming codes. 



Generalization to the qudit case: phase matrices In the following we generalize 
the approach described in [I] to pair-interactions between higher-dimensional systems, i.e., 
qudits. It will be useful to express a general pair- interaction Hamiltonian with respect to 
a so-called nice error basis and to use the matrices from such basis as control operations. 
First, we recall the definition of nice error bases |15 | I24 | 125]. 
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Definition 3 (Nice error basis) Let G be a group of order d 2 with identity element e. 
A nice error basis on C d is a set £ = {U g G C dxd \ g G G} of unitary matrices, which are 
labeled by the elements of G, such that (i) U e is the identity matrix, (ii) trU g = d 5 g ^ e for 
all g G G, and (in) U g Uh = a(g,h)U g h for all g,h G G. The factor system a(g,h) is a 
function from G x G to the set C x := C \ {0}. 

Condition (ii) shows that the matrices U g are pairwise orthogonal with respect to the trace 
inner product. The group G is called index group since its group elements index the elements 
of the nice error basis £. 

Example 4 Let d G N and let uj = exp(2iri/d) denote a primitive d-th root of unity. Next, 
we define operators S := X^=ol^)(^ + 1|> w ^ iere ^ e indices are reduced modulo d, and 
T := Efc=o' jfc | /c )( /c l- Then the setS d ■= {5*2* : i = 0,...,d- l,j = 0,...,d- 1} is a 
nice error basis on C d (see, e.g., )2b^ ). This shows the existence of nice error bases for 
any dimension d G N. In this case the index group is the abelian group G = 7Ld x Z^. The 
identity ST = toTS is readily verified. This shows that the corresponding factor system a 
is given by a((i, j), (k, £)) = u)~i k , for all (i, j), (k, £) G G. 

Next, we describe a particular way of representing a general pair-interaction Hamiltonian 
acting on n qudits. Let £ = {U g : g G G} be a nice error basis. Since the matrices U g form 
a basis of <C dxd a general pair-interaction Hamiltonian on n coupled qudits may be written 
as 

k<i h,h'^e k h^e 

where the coefficients Jki-hh' £ C and G C are chosen such that H is a traceless Hermi- 
tian matrix. In the following we consider decoupling schemes with time-slots of equal length 
and with elements of a nice error basis having abelian index groups as control operations. 
Then, in each time-slot, each Uh in eq. Ijlljl acquires a phase factor that is controlled by 
the applied local unitaries of the nice error basis. We define x(g, h) to be the phase factor 
that Uh acquires when it is conjugated by U g , i. e. x(<5S h) is defined via the relation 

ulu h u g = x(g,h)u h . (12) 

Hence, the bilinear term Uh®Uh> acquires the phase factor x(d, h)x(g' ', h') if it is conjugated 
by Ug®U g i. The fact that conjugation by U g of Uh merely introduces a phase factor follows 
from the index group being abelian. We need to characterize the corresponding d 2 x d 2 
matrix 

X:=(x(9,h)) 9iheG (13) 

with entries x{di h). 

Lemma 5 Let £ := {U g \ g G G} be a nice error basis with an abelian index group G. Then 
the matrix X = (x(g 5 h)) g h&G * s ^ e character table of the group G. 
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We give a proof in Appendix El This lemma shows that X is a character table of G 
and thereby provides a natural generalization of Table El to non-qubit systems. Given a 
decoupling scheme which contains operators from a nice error basis £ = {U g : g £ G} we 
obtain a collection of phase matrices Ph for each h £ G in the following way: the (fc,j)th 
entry of Ph — denoted by Ph;k,j — is given by the phase factor of Uh that is acquired in the 
jth time-slot by conjugating the k qudit by the respective element of the scheme. 

Criteria for decoupling in term of phase matrices Let Ph be phase matrices. Then 
they define a decoupling scheme if and only if 

N 

Yl = (14) 

3=1 

for all k and for all h ^ e, and 

N 

Ys P h;kjPh';ej=0 (15) 
i=i 

for all k < £ and for all h, h! ^ e. The condition given in eq. (|14j) ensures that all local 
terms are removed. Eq. (|15fl implies that all bilinear terms are removed. The question how 
to construct collections of phase matrices which satisfy eqs. (j!4j) and (|15|) is addressed in 
Section 15.11 and a solution is presented if all nodes have a dimension which is a power of a 
prime. 

4.3 Decoupling schemes from orthogonal arrays 

Orthogonal arrays have been applied in the design of experiments to plan statistical data 
collections systematically. The books HH1 HZ) provide good introductions to the topic. 
In the following we recall the definition of orthogonal arrays (or OAs for short) and show 
how to use them for the decoupling problem. M. Stollsteimer and G. Mahler have first used 
OAs for the construction of decoupling schemes and selective coupling schemes for pair- 
interaction Hamiltonians acting on qubits. This method was generalized to pair- interaction 
Hamiltonians acting on qudit 5 j with the helps of unitary error bases. Even more generally, 
we show in this section how to use orthogonal arrays for constructing decoupling schemes 
for t- local Hamiltonians acting on qudits with t > 2. 

Definition 6 (Orthogonal array of strength t) Let A be a finite set and let n,N G N. 

An n x N array M with entries from A is an orthogonal array with \A\ levels, strength 
t, and index A if and only if every t x N sub-array of M contains each possible t-tuple of 
elements in A precisely A times as a column. We use the notation OA\(N, n, s, t) to denote 
a corresponding orthogonal array. If X, s, and t are understood we also use the shorthand 
notation OA(N,n). 

In statistics for the various parameters of an orthogonal array (OA) some traditional 
terminology is used. In the context of decoupling and simulation of Hamiltonians different 
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Parameter 


Design of Experiments 


Decoupling Schemes for qudit systems 


n 


factors 


nodes (qudits) 


N 


runs 


time-slots 


A 


levels 


elements of an operator basis 


d 


number of levels 


(dimension of the nodes) 2 


t 


strength 


locality 


A 


index 





Table 1: Dictionary between notions used in the theory of design of experiments and the 
theory of qudit systems to describe the parameters of an orthogonal array OA\(N,n,d,t) 
over alphabet A. 



names for these parameters are used. We provide a dictionary between the different lan- 
guages in table ^ Note that as a convention we write OAs asnxJV matrices, whereas 
most authors in design theory prefer to write the matrices as N x n matrices. Besides 
typographic reasons we found the presentation using n x N matrices useful to establish the 
correspondence with pulse sequences in NMR which are typically read from left to right like 
a musical score |10| . 

An important special case arises if the strength t is two. This means that each pair of 
elements of A occurs A times in the list ((akj,aej) I J ' = 1> • • • N) for 1 < k < I < n. Most 
of the known constructions actually yield arrays of strength two (22,- For many physical 
systems it will be sufficient to study arrays of small strength since the strength relates to the 
degree of the interactions, i.e., for pair-interaction Hamiltonians it is sufficient to consider 
arrays of strength t = 2. 

Example 7 As an example of small size we give an orthogonal array with parameters 
OA{lQ, 5, 4, 2). This means that we have 16 runs/time-slots, 5 factors/qubits, 4 different 
symbols/pulses, and (pair-interaction) strength two. 

/1111222233334444\ 
1234432112344321 
1342134231243124 
1342421342131342 

\1234214343213412/ 

It is straightforward to check that indeed every pair of rows contains all the 16 possible 
pairs of symbols precisely once. This array was obtained from a linear error-correcting code 
over the finite field F4. We will explore this construction in more detail in Theorem I1UI 

Next, we generalize the ideas in |S] and describe how to achieve decoupling of t- 
local Hamiltonian acting on qudits. The basic idea is to use an orthogonal array M with 
parameters OA(N,n,d 2 ,t) over an alphabet A of size d 2 . Here d denotes the dimension of 
the qudits. The elements of A are identified with the operators Ui,... ,U d 2 of a unitary 
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error basis [3]. The orthogonal array is an n x N matrix M = (rrikj) and determines which 
elements of the unitary error basis are used on the qudits in the time-slots as follows: 

1 - 

n • • • ® U U H ( U ^ ® • • • ® (is) 

3=1 

The following theorem shows that the resulting time evolution is that of the zero Hamilto- 
nian which means that indeed the decoupling conditions given in Definition ^ are satisfied. 

Theorem 8 (Decoupling with OAs) Any orthogonal array OA(N,n,d 2 ,t) over an al- 
phabet of size d 2 can be used to decouple n qudits which are governed by a t-local Hamiltonian 
within N time-slots. 

Proof. First, note that for any t-local Hamiltonian H of a system consisting of t qudits the 
following operations define a decoupling scheme jSJ: Let Ux, U2, ■ ■ ■ , U d 2 denote the elements 
of some unitary error basis for C°\ Since the tensor products of all possible pairs of these 
elements form a vector space basis of the linear maps acting on (C d )^ we obtain that [12 IS] 

1 S 

— ]T (tft ®...®ul)H(U h ®...®U h ) = (17) 

ii, — ,*t=l 

for all (traceless) Hamiltonians acting on (C d )®*. Recall that we assume without loss of 
generality that H is traceless. Let B := {a\, 02, . . . , cr d 2_ 1 } be a basis for the vector space 
of traceless matrices of size d x d. Recall that a general t-local Hamiltonian acting on n 
qudits can be written as 

t d 2 -i 

H = E E E J i-. .,o-;!- : ...o-:t- (is) 

s=l (ki,...,k 3 ) a 1 ,...,a s =l 

where the second sum runs over all s-tuples with different entries from {l,...,n} and 
J(k 1 ,...,k B );a 1 ,...,a s £ C. Now, we pick any s-subset {ki, . . . , k s } C {1, . . . , n} of the nodes and 
denote by Ck u ,..,k s the coupling between these nodes. We define Cfc li ... j fc s to be the coupling 
among the qudits ki, ■ ■ ■ , k s , i.e., 

C kl ,...,k. ■= E ■ 7 (fci,...,fc.)iai,..,a. ® • ■ ■ ® CT £ S) • 

We define ___^ s to be the corresponding operator acting on (C d )® s . Formally, we say that 
Cfci,...,fc s is obtained by embedding into (C^)®" - according to the tuple (fei, . . . , fc s ). 

For any operator X acting on (<£ d )® s we denote the embedding into (C d )® n according to 
the tuple (k u ...,k s ) by X^'-' ks \ 

The idea of the proof is to reduce the problem to eq. (|17|) by using the local structure 
of the Hamiltonian. Since M is an OA(N, n, d 2 ,t) all elements of {1,2,..., d 2 } 3 for s < t 
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appear equally often in the list (jn klt j, . . . , mk a ,j) where j = 1, . . . , N. Therefore, the average 
Hamiltonian corresponding to the coupling among the qudits fci, . . . , k s is evaluated as 
follows: 



1 N 
i=i 
1 N 



■®Ul k JC kl ^ ks (U mki j 



u„ 



3=1 



(ki,...,k 3 ) 



d 2 



d 2t 



u is ) 



ii,...,i s =l 



(ki,...,k s ) 



0. 



The equality between the second last and last line is due to eq. (fTTj) . 



□ 



5 Equivalence of the Constructions 

We show that the methods based on phase matrices and orthogonal arrays of strength two 
lead to the same class of decoupling schemes if we use elements of a nice error basis with 
an abelian index group as control operations. More precisely, we prove that the decoupling 
conditions given in eqs. (|14|) and (|15j) are equivalent to the condition that the decoupling 
matrix is an orthogonal array of strength two. 

5.1 Phase matrices from orthogonal arrays 

We show that a decoupling scheme based on an orthogonal arrays of strength two gives rise 
to phase matrices satisfying the decoupling criteria eqs. (fTl|) and (fH)|) for pair- interaction 
(two-local) Hamiltonians. 

Theorem 9 Let G be a finite abelian group of exponent e(G) and let £ = {U g \ g £ G} be a 
nice error basis for C d with index group G. Then a decoupling scheme for n qudits governed 
by a pair-interaction Hamiltonian can be constructed. The scheme uses N time-slots and 
is constructed from an orthogonal array OA(N, n, d 2 , 2) over the alphabet A = {1, 2, . . . d 2 }. 
Furthermore, the scheme gives rise to phase matrices [PhjheG of size n x N with entries 
in {1, lj, . . . , oj e ^ G ^ 1 } which satisfy the orthogonality conditions and are compatible with 
respect to taking Schur products. 

Proof. We denote the orthogonal array by M = [wiftj], where k = 1,... ,n and j = 
1, . . . , N. Next, we fix an ordering g\, ... , g d 2 of the elements of G and assume that g± = e 
is the identity. Moreover, we identify the operators of £ with the elements of A according 
to 1 i — > Q\ , 2 g 2 , . . . ,d 2 i — > g d 2. Note that conjugating Uh £ £ by U g £ £ results in a 
phase factor that is the (g, h)th entry of the character table of G. 
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Starting from the given orthogonal array we construct d 2 phase matrices P gi , . . . , Pg d2 
as follows. Pick row number k of the OA and replace each symbol according to a i— > v 9a 
for a = l,...dr where v ga denotes the g a th row of the character table of G (cf. Lemma 0). 
By assigning P g% for i = 1, . . . , d 2 the first, second, and d 2 th components of each entry, we 
define the rows number k of the d 2 phase matrices. In other words, the entry (k,j) of the 
matrix P^, where h £ {gi, . . . ,5^2}, is given by Ph;k,j = x( m k,j, h). Note that the matrix 
P gi is the all-ones matrix of size n x N. 

While the condition P g oP h = P gh is automatically guaranteed since the characters form 
a group, we have to show that the resulting vectors are pairwise orthogonal. In order to do 
so we pick two rows k and £ of the original orthogonal array. We may assume that the two 
rows have the following form (or else we apply a suitable permutation of the columns) 



1 1 
1 2 



1 

d 2 



2 2 
1 2 



2 

d 2 



d 2 



d 2 
2 



d 2 
d 2 



(19) 



A times 



since all pairs appear equally often (A times) in the OA. Let A = (+ + •••+) be the vector 
of length A = N/d 4 containing only the entry +1. Furthermore, for each h £ {g±, . . . ,<7d 2 } 
define a vector Wh G C^ 2 as follows. We define Wh '■= [x(<7i, h), . . . , x{9d 2 i i- e -> w h is the 
hth column of the character table of G (cf. Lemma ISJ). By substituting the entries of rows k 
and £ in the form of eq. (|19j) by the corresponding sign assignments in Table © we obtain 
as the kth rows of P, 
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, Pg d2 the vectors 



A (8) w gi <S> (H h) , A (8) w g2 (g) (H h) , 

and for the £th rows of P Sl , . . . , P 9d2 the following vectors: 



A <g> w, 



9d* 



(+•••+) 



A ® (H h) <8> , A (8) (H h) (8 tu. 



92 ' 



A (8 (+•••+) <8 w 



9 d z 



where (+•••+) is a vector of length d 2 . Whenever gi, gj are not both equal to the identity 
gi all these vectors are orthogonal to each other since the columns of the character table 
are orthogonal. This shows that all rows of the matrices P g2 , . . . , P 9d2 are orthogonal and 
the matrices satisfy the Schur condition P g o P h = P gh . □ 



Phase matrices from Hamming codes In the following we show how the known con- 
structions of sign matrices can be reproduced with well known families of orthogonal arrays. 
What is more, we show that the class of orthogonal arrays used for this construction are 
based on Hamming codes |28M29| . They can be used to construct phase matrices for higher- 
dimensional systems in case the dimension d of the qudits is a power of a prime. 

First, we briefly recall some basic facts about error-correcting codes |281 129j since they 
will feature in the subsequent constructions of orthogonal arrays. A linear code over the 
finite field ¥ q is a fc-dimensional subspace of the vector space F™. The metric on the space 
F^ is called the Hamming weight. For x = (xi,...,x n ) £ F^ we have that wt(x) := 
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\{i G {1, ...,n} : Xi ^ 0}|. The minimum distance of a linear code C is defined by 
d = dmin := min{wt(c) : c G C, c 7^ 0}, where denotes the zero vector. As a shorthand 
we often abbreviate this situation by saying that C is an [n, k, d] q code. We need one more 
definition which is the dual code C 1 - of C defined by C 1 - := {x G F™ : x-y = Ofor ally G C}. 

The following theorem |27l Theorem 4.6] establishes a connection between orthogonal 
arrays and error-correcting codes. In fact this is one of the most prolific constructions for 
OAs known. 

Theorem 10 (OAs from linear codes) Let C be a linear [n, k, d] q code over ¥ q . Letd 1 
be the minimum distance of the dual code C ± . Arrange the code words ofC into the columns 

k 

of a matrix A G ¥ q xq . Then A is an OA(q k , n, q, d 1 - — 1). 

For the case of a network consisting of n qubits which are governed by a pair-interaction 
Hamiltonian we can construct decoupling schemes using N pulses from an OA(N, n,4, 2). 
Hence, in order to apply Theorem 1101 we have to find a code C of linear codes over F4 for 
which the parameters are [n, k, d] and for which the minimum distance d 1 - of the dual code 
is at least 3. 

Let q be a prime power and let m G N. Then the Hamming code H qtUl of length 
n = {q m — l)/(q — 1) is a single-error correcting code over the field ¥ q with parameters 
[n,n — m,3]q. The dual code H q ^ m of the Hamming code [n,n — m,3] g has parameters 
[n,m,q m ~ 1 ]. By specializing q = 4 and by using Theorem 1101 for H^ m we therefore obtain 
orthogonal arrays with parameters OA(N, n, 4, 2), where n = (4 m — l)/3 and = 4 m for 
any choice of m G N. The alphabet set is in this case the finite field F4 of four elements. 

The procedure to obtain a decoupling scheme for a network of n qubits, where n is an 
arbitrary natural number, i. e., not necessarily of the form n = (4 m — l)/3 is as follows: first 
let m G N be the unique integer such that n < 3 < 4n. Then construct the orthogonal 
array with parameters OA(4 m , (4 m — l)/3,4,2). The columns of this OA are code words 

I (4 m — 1)/3 

of Hf m C F4 . We can now obtain a triple of sign matrices S x , S y , and S z by using 

the substitution rules in Theorem |SJ This leads to the same sign matrices as the ones 
constructed in [I] by a direct construction and in jS] using spreads in the geometry F| m . 

In case the dimension is an arbitrary power of a prime d = p r , we use the Hamming 
code [n, n — m, 3]^2 to obtain an OA(N, n, d 2 , 2), where n = (q m — 1) /(q — 1) and N = q m . 
By Theorem El from this orthogonal array a collection of phase matrices can be constructed 
which satisfy the orthogonality conditions and are compatible with respect to taking Schur 
products. 

5.2 Orthogonal arrays from phase matrices 

In this section we provide a converse to the previous section by showing that orthogonal 
arrays of strength two can be constructed from phase matrices satisfying the decoupling 
conditions in eqs. ()14|) and H15|) . To do this we need the following lemma which gives a 
criterion in terms of group characters to decide whether an element of the group ring is an 
equally-weighted sum of all group elements. This allows to check whether a matrix is an 
orthogonal array. 
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Lemma 11 Let G be an abelian group of order \G\ . Denote by xi,X2, ■■■ , X\G\ a ^ irreducible 
characters of G, where xi is the trivial character (i. e. Xi(h) = 1 for all h E G). Let v 
be an arbitrary element of the group ring C[G], i. e., v is a formal sum of (weighted) group 
elements 

v:=^2fi g g, |J 9 eC. (20) 

g eG 

If Xi{ v ) = for all i = 2,...,\G\ then we have v = -^^2 g£ Q g , where fi := Xi{ v ) = 

A proof of this lemma is given in Appendix FBI We are now ready to state the main 
result of this section. 

Theorem 12 Let G be a finite abelian group and let \P g \ g ^G be a collection of phase 
matrices of size n x N which are compatible with respect to taking Schur products and 
satisfy the orthogonality relations. Then these phase matrices define an orthogonal array 
OA(N,n,\G\,2). 

Proof. For fixed k = 1, . . . , n and j = 1, . . . ,N each vector Vkj ■= [Ph] k, j]h^G is a row the 
character table of G. Therefore, it determines uniquely g E G such that the entries of v^j 
are Xg(h), i. e., the values of the irreducible character corresponding to g applied to h G G. 
We denote the so defined group element by g^j- Let M = (gkj) be the n x N matrix with 
entries gkj. We will show that M is an orthogonal array OA(N,n, |G|,2). 

Pick any two rows (gkj) and (gij) of M. We define an element of the group ring C[GxG] 
as the formal sum 

N 

ru ■= /A9kj,9lj) ■ 

3=1 

To abbreviate the notation we denote by Xg,g' ^ ne irreducible character of G x G corre- 
sponding to the element (g,g'). The decoupling conditions given in eqs. (|14|) and l|15j) are 
equivalent to 

Xg,g'{rkl) = 

for all (g,g') ^ (e, e). By Lemma 1111 this is equivalent to the case that all elements of 
G x G appear equally often in the sum r^. This shows that M is an orthogonal array 
OA(N, n, \G\ , 2) of strength t = 2 over G. □ 

Theorems El and El shows that phase matrices can be used to define an orthogonal 
array of strength two and vice versa. Based on the above lemma we give an alternative 
characterization of orthogonal arrays of arbitrary strength which is a generalization of 27, 
Theorem 3.30]. This theorem implies that the entries of the array can be replaced by 
complex numbers such that the resulting matrix is orthogonal with respect to the usual 
inner product for strength t = 2. Recall that for elements of the Cartesian product v E G n 
the Hamming weight wt(v) is defined by the number of components which are different from 
the identity. 
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Theorem 13 (Conditions for a matrix to be an OA) Let G be a finite abelian group 
and let A be a matrix of size n x N with entries from G. Then A is an orthogonal array 
OA(N,n, \G\,t) if and only if 

N n 

EIIx(AX> = (21) 

3=1 i=l 

holds for all nontrivial characters \ £ Irr(G) and for all v £ G n of Hamming weight 
•wt{v) < t. 

Proof. Suppose A is an orthogonal array OA(N, n, \ G\,t) and let v £ G n be a fixed element 
of weight wt(v) < t. Let be a non-trivial character of the t-fold direct product G t . Denote 
by vt £ G t the vector containing the components of v which are different from the identity. 
We obtain that ^2 g&G t 'O(gvt) = 2~^ g eG* ^(d) = anc ^ * ne statement follows from the fact 
that the characters of G t are given by products of t characters of G. 

Conversely, assume that we are given a matrix A such that cq. Ij21j) holds for all non- 
trivial characters and all v £ G n of weight wt(v) < t. In particular, this means that eq. 1)2 1[) 
is satisfied for the vector w e all components of which are equal to the identity e £ G. Again, 
fix a t-subset T of the rows of A. We have that for all non-trivial characters i? of G l the 
identity Y^,j=i ${9jW e ) = holds, where the elements gj £ G* are obtained by selecting the 
jth column of A, followed by selecting the components corresponding to T, and finally to 
consider the element as being an element of G . Now, we can apply Lemma ITT1 to obtain 
that the list [gj : j = 1, . . . , N] has to contain all elements G t and that each element has to 
occur the same number of times. This shows that A is an OA(N,n, \G\,t). □ 

5.3 A new characterization of orthogonal arrays 

In Section 14.31 we have seen that orthogonal arrays of strength t can be used to construct 
decoupling schemes for t-local Hamiltonians. In order to establish a converse result we need 
some additional conditions on the class of schemes considered: (i) the schemes have to be 
regular (see Definition ^) and moreover we will assume that (ii) the pulses are actually taken 
from a fixed set of unitaries which in addition will be assumed to form a unitary error basis. 
We begin by stating some standard concepts from quantum information theory which will 
be used in the proof. Recall that the Shannon entropy is defined by the equation 

M 

H{pi ,...,p M ) = -J2 p i log 2 Pi ■ 

3=1 

Shannon entropy measures the disorder of probability distributions. If pj = 1 for some j, 
then the entropy is zero. The entropy takes its maximum value log 2 M for the uniform 
distribution. The notion of entropy extends to density operators, and is usually called von 
Neumann entropy. Let p be an arbitrary density operator on C d . Then the spectral decom- 
position p = Ylj=i^j\^j}(^j\ is such that the eigenvalues Ai,...,Aa/ form a probability 
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distribution and the eigenvectors • • • , \ ^m) form an orthogonal basis of C M . The von 
Neumann entropy S(p) of p is defined by the equation 



M 

S{p) = -^Ajlog 2 Aj 



The von Neumann entropy takes its minimal value on pure states, i. e., for p = |^'}(^ r |, and 
its maximal value log 2 M for the maximally mixed state p = 1/M. Let U±, . . . , Un £ (£ MxM 
be arbitrary unitary matrices, pi, . . . ,pn a probability distribution, and \^>) a state of C M . 
We have the following inequality (see [3U], p. 518) 

S (X>^]l*>W'j <H(p u ..., PN )<log 2 N. (22) 

The following theorem shows that if a regular decoupling scheme can switch off an arbitrary 
t-local Hamiltonian, then the tensor products of the operations performed on an arbitrary 
t-tuple of qudits must form a unitary error basis for this subsystem. 

Theorem 14 (Equivalence of decoupling schemes and OAs) 

Let D be a regular decoupling scheme that uses elements of a unitary error basis £ := 
{U\, . . . , Utf} as control operations, acts on n qudits and consists of N time-slots. Denote 
by U mij (g) U m2j (g) • • • (g> U mnj the local operation that is performed on the qudits in time- 
slot j = 1,...,N, i.e., the indices m^j E {1, . . . , d 2 } determine which elements of £ are 
applied to the qudits in the time-slots. The scheme D can be used to decouple any t-body 
Hamiltonian if and only if the matrix M = (m^j) where k = 1, . . . , n and j = 1, . . . , N is 
an orthogonal array OA(N, n, d 2 ,t) of strength t. 



Proof. Consider a fixed t-tuple (ki, . . . ,kt) with different entries from {1, . . . , n}. Let 

H k ,.., kt :=E E E ^wWS' • • • , 

8=1 {£i,...,£ s ) a i- a * 

where (t\, . . . , £ s ) runs over s-tuples with different entries from {k±, . . . , kt}- We say that 
the operator ^ is the restriction of the t-body Hamiltonian H to the qudits ki,...,kt- 

We denote by Hk 1: ... : k t the corresponding operator acting on (C d )®*. 

Note that for every traceless Hermitian operator X acting on (C d )®* there is a i-local 
Hamiltonian H such that its restriction Hfa ... & ( to the qudits k\,...,kt is given by the 
embedding x( kl >-> ki > of X to (C d ) 0n according to the tuple (ki,...,kt). Let T D be the 
operator 

N 

T D : ^2pj(U miJ ® U m2j ® ■ ■ ■ ® U mn] ) ] H{U mi] U m2j ® • • • ® U mnj ) 

i=l 



19 



Define the weight i% of each tuple (i±, . . . , it) G {1, . . . , d 2 } f to the sum of all pj's with 
( m k!,j, ■ ■ ■ j mksj) = (h, ■ ■ ■ , k)- Now suppose that Td{H) = for all t-body Hamiltonians. 
Consequently, we have that T£>{Hki ... k t ) = f° r an restrictions to i-tuples. But this implies 
that the weights for all A 1 must be equal. This is seen as follows: the equality 



T D (H, 



k\...k t 



w ii...h( U h ® 

i\,...,it=l 

Q(fc 1 ,...,fct) _ Q 



(fci,...,fc t ) 



shows that the operation defined by the sum above is a unitary depolarizer for (C )®*, i. e., 
( / 2 



X] ^1.-^(^1 ® •••®?7iJ t ^(C/ il ® ■ ■ ■ ® U it ) 



ii,...,it=l 



tr{X). 
dP 



for all operators X acting on (C d )®'. Now, we show that the weights must be all equal. Let 
• • • , be an orthonormal basis of (C rf )®*. We define a special state in the bipartite 
system (C d )®* (g> (C d )®* together with its corresponding density operator 



^= rft S l*r>(*.|® |#r>(*< 



r=l 



r,s=l 



We use the fact that Tjj is a unitary depolarizer to show that show that all weights are 
equal 

l®T D {p) = Y w il ... lt (l d t®Ui 1 ®---®Ui t ) J! p(l d t<g)U ll <g)---®Ui t ) 
h...i t eA t 

= M E l*r><*.|® 



r,s=l 



X] (tfii ® • • • C/ it) t l*r)(* S | (^i 

ii...i t e^.* 

1 d 

-£|¥ r >(¥ r |®:Wd* 



E4 



r=l 



2/ 



= l rf t /d* <g> \ A t jd 1 = l d 2t jd 

It follows from the above equation that we need at least d 2t different unitaries are necessary 
since the rank of each pure state (1 ® ?7^)| X I , )(^ / |(1 (g> U) is one and since they have to sum 
up to d 2t (the rank of the maximally mixed state) . 

Since we use exactly d 2t different unitaries (tensor products of elements of the unitary 
error basis £) as control operations all weights u^...^ must be equal due to the inequal- 
ity (|22|). Now together with the fact that for regular schemes all time-slots have equal 
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length we conclude that (g> • • • <£>Ui t must appear with the same multiplicity. Therefore, 
by considering all t-tuples k\ , . . . , kt oft qudits we see that the decoupling D scheme must 
correspond to an orthogonal array OA(N, n) with d 2 levels and strength t. □ 



6 Conclusions 

We have shown the equivalence between two constructions for decoupling schemes and 
selective coupling schemes in networks of qubits. One construction is based on triples of sign 
matrices which are closed under taking entry- wise products, while the other construction is 
based on orthogonal arrays of strength two over an alphabet of size four. 

The construction using orthogonal arrays can be generalized to systems where the nodes 
have higher dimensions. Also the case where the system Hamiltonian has higher couplings 
can be dealt with by using orthogonal arrays: the coupling order directly translates into 
the strength of the orthogonal array. A special case arises when the local pulses which are 
applied in each time-slot are actually elements of a nice error basis for an abelian group. 
We have shown that in case of equidistant interval lengths (after refinement) this leads to a 
class of schemes which are equivalent to orthogonal arrays. In addition we have presented 
a construction of schemes for decoupling and selective coupling which can be constructed 
by using Hamming codes. 

Moreover, we have shown that the construction of this particular class of decoupling and 
coupling schemes can be reduced to questions about the existence of these combinatorial 
arrays. While several constructions for orthogonal arrays are known, there still remain 
some open problems such as the case where the dimensions of the nodes could be different. 
Another important problem is to devise schemes for a situation where the given Hamiltonian 
is of a particular form, i.e., where not all interactions are present or can be assumed to 
be very weak for a large number of pairs. In this case a combination of graph theoretical 
techniques and the methods described in this paper can be developed. 
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A Appendix: Characters of abelian groups 



In this appendix we collect some basic facts of the representation theory of finite groups 
which are needed in the paper. Recall that GL(n, C) denotes the group of invertible n x n 
matrices with entries in C. Let G be a finite group. A representation of G over C (see also 
[^E21) is a homomorphism p from G to GL(n, C), for some n G N. The degree of p is given 
by n. The representation p is called irreducible if there are no invariant subspaces under 
the action of the matrices {p(g)} g ^G apart from the trivial subspaces {0} and C n . With 
each n x n matrix p(g) we associate the complex number given by the trace of p(g), and 
call this number x(d)- The function x '■ G — > C is called the character of the representation 
p. A character is called irreducible if the corresponding representation is irreducible. We 
need the following theorem on abelian groups [A2\ Chap. V. §6]. 

Theorem 15 (Characters of abelian groups) Let G be a finite abelian group of order 
\G\. Then every irreducible representation p of G has degree 1, i. e., we have that p : G — > 
C x is a homomorphism which maps G to scalars. Furthermore, the number of different 
irreducible representations (irreducible characters) of G is given by \G\ and the characters 
form a group G = Hom(G, C x ) under pointwise multiplication. Hence, we have that 

XX(h) = x(h)x(h) 

for all irreducible characters x> X ana " h £ G. Moreover, the character group G is isomorphic 
to G. Thus, we can label the characters of G by the elements of G using an isomorphism 
which maps ft i— » Xh for all h £ G. 



B Appendix: Proof of Lemmas and 1111 

Lemma |SJ Let £ := {U g \ g £ G} be a nice error basis with an abelian index group G. 
Then the matrix X = (x(d, h)) g h£G is the character table of the group G. 

Proof. Let a be the factor system corresponding to the nice error basis £ with abelian index 
group G. We prove that X is a character table by showing that the rows of X form a group 
under pointwise multiplication that is isomorphic to G (see Theorem 1151 in Appendix). We 
first show that 

a(h,g) 

X(g,h)= (23) 
a{g,h) 

We have that 

UgU h = a(g,h)U gh , (24) 
U h Ug = a(h, g)U hg = a(h, g)Ug h . (25) 

By multiplying eq. (|25|) by Ug from the left and using eq. (|24j) we obtain 

U\U h Ug = a{h,g)U\U gh 



24 



= ^4utu 3 u h 

a(g,h) g 
a(g,h) 

We now prove that the rows of X form a group under pointwise multiplication that is iso- 
morphic to G. Let g, g be arbitrary elements of G. Note that we have a g) a (g , g) = 1 
(otherwise the matrix U g -iU g = a(g , g)U^-x g would not be unitary). The group property 
is verified by 

x(g,h) x (g-\h)u h = ulu^UhU-g-iUg 

= «(5" 1 ; 5)a(5" 1 ,5 f ) U lg-l U hUgg-l 

gg 1 ™ 99 

for all h E G. 

The rows of X form a group that is isomorphic to G (and not only to a proper subgroup 
of G) since there is a bijection between the rows of X and the elements of G. This is seen 
as follows. Assume that there are g ^ g such that x{9i^) = xid^h) for all h E G. This 
is equivalent to V\\J h Ug = U ] g U h Ug. Set U = UgU\. Then we have UM = MU for all 
M E C. dxd since the matrices form a basis of C dxd . Therefore U must be a multiple 
of the identity matrix. Due to the properties of a nice error basis this is only possible for 
g = g. This proves that there is a bijection between the group elements of G and the rows 
of X. □ 



Lemma II U Let G be an abelian group of order \G\. Denote by Xli X2> ■ ■ ■ > X\G\ °ll irre- 
ducible characters of G, where xi is the trivial character (i. e. Xi{h) = 1 for all h E G). 
Let v be an arbitrary element of the group ring C[G\, i. e., v is a formal sum of (weighted) 
group elements 

v:=^^ g g, /j, g eC. (26) 
g&G 

If Xi( v ) = for all i = 2,...,\G\ then we have v = r^rYlgeG 9 > wnere I 1 := Xi{ v ) = 

Proof. Let G := {g±, . . . ,g\G\} be an arbitrary ordering of the group elements, where g\ is 
the identity element of G. Denote by X the (normalized) character table of G, i.e., 

Xij'=\G\- 1,2 Xi{9j) (27) 

for i, j = 1, . . . , \G\. Recall that the (normalized) character table X is a unitary matrix and 
has the following form |32l 13 1 j 

/ 1 ... 1 



x 1 



|G|V2 



: * 

V 1 



(28) 
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The conditions given in the lemma can now be expressed as 

/' 


V o / 



lei 1 / 2 * 

V M|G| / 

Multiplying by the inverse X~ l we obtain 



/ A*i \ / A* \ 

^2 



due to the special form in eq. (|28|). This show that all coefficients /U 9 in eq. 1)20(1 are equal 
to □ 
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